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I  INTRODUCTION 


A  critical  problem  in  a  missile  defense  system  is  that  of  predicting 
the  trajectory  and  impact  point  for  a  ballistic  reentry  vehicle.  This 
memorandum  will  describe  methods  of  prediction  as  well  as  numerical  results 
for  several  representative  examples.  There  are  several  reports  describ¬ 
ing  the  estimation  of  the  states  of  a  ballistic  missile;  the  ballistic 
trajectory  and  impact  point  will  be  predicted  by  using  these  estimated 
values. 

There  are  several  important  points  to  be  considered.  First  is  the 
choice  of  the  coordinate  system  to  be  employed.  Cither  a  radar  coordinate 
system  or  a  rectangular  coordinate  system  centered  at  the  radar  site  can 
be  used  for  the  problems  being  studied.  Second  is  the  treatment  of 
physical  parameters  in  the  equation  of  motion,  such  as  the  ballistic  co¬ 
efficient  of  the  reentry  vehicle  and  the  eccentricity  and  rotation  of  the 
earth.  Since  the  time  required  for  computation  may  become  significantly 
large,  it  is  also  very  desirable  to  find  a  closed-form  solution  of  the 
equation  of  motion,  which  is  a  rather  complex  nonlinear  differential 
equation.  The  important  point  here  is  how  much  the  accuracy  of  the  solu¬ 
tion  is  degraded  in  obtaining  a  closed-form  solution.  Third  is  the  propa¬ 
gation  of  initial  errors  to  the  final  values  in  prediction;  this  is  val¬ 
uable  in  order  to  trade  off  the  magnitude  of  errors  and  the  computation 
time  in  estimation.  For  this  purpose,  the  sensitivity  of  the  initial 
values  to  impact  points  is  briefly  investigated. 
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in  which 


[x,y,s,x,y,z]  *  target  state 
t  rth  rotation  rate 

Product  of  gravitational  constant  and  mass 
of  earth 

Geodetic  latitude  of  radar 

Geocentric  latitude  of  radar 

Atmospheric  mass  density  at  target  position 

Ballistic  coefficient  (t.e.,  weight -to-drag  ratio) 

Earth  radius  to  radar  site 

Magnitude  of  position  vector  from  earth  center 
to  target 

Velocity  magnitude  of  target 
Gravitational  acceleration 

f,  and  are  the  initial  values  of  r.  V  and 


4 


f  =  Ki-riMit  r  it- i  l  y  o!'  re  f  “tence  ellipsoidal  earth 
0.0066945) 

J  Dimensionless  constant  -  1.624  x  10  ^ 

■t  ~  Geocentric  latitude  ol'  target 
a  -  Kcpiutort.il  radius  of  earth  (20,926,743  ft). 

'Hi e  linear  coeffic  ient  matrix  A  is  a  6  *  6  matrix  whose  elements  are 
constant  except  for  three  .  leirents  containing  the  atmospheric  mass  density 
The  vector  B  is  constant  and  the  vector  C  is  a  nonlinear  term  that 
•a  1  1  1  be  negligible  if  T,  the  magnitude  of  the  target  velocity,  and  r,  the 
magnitude  of  the  position  vector  from  the  earth  to  the  target,  do  not 
change  significantly.  The  last  term  D  contains  the  elements  describing 
the  influence  of  the  eccentricity  of  the  reference  ellipsoidal  earth. 

Hence,  if  eccentricity  *•  is  considered  to  be  lero,  then  the  tern:  D  vanishes. 

One  of  the  objectives  of  this  report  is  to  investigate  simplifications 
of  the  differential  equation  described  above.  If  the  terms  C  and  Q  are 
negligible,  the  differential  equation  will  become  A  -  4A  ♦  8.  It  is  true 
that  a  linear  differential  equation  with  time  varying  coefficients  is  no 
better  for  finding  an  analytical  solution  than  a  nonlinear  differential 
equation.  However,  if  the  time  varying  coefficients  are  approximated  as 
constant  for  a  certain  time  interval,  then  piecewise  closed-form  solutions 
can  be  obtained. 

Based  on  physical  considerations,  it  is  helpful  for  the  purpose  ol 
the  following  discussion  to  divide  the  atmosphere  into  two  regions.  One 
called  eiouf eovphere,  defined  as  the  space  above  an  altitude  of 
300.000  ft;  the  other  i»  called  endooMosphere,  defined  as  the  space  below 
an  altitude  of  .300,000  ft. 
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Ill  ENDOATMOSPHERIC  PREDICTION 


For  ballistic  missiles  at  an  altitude  of  300,000  ft  with  a  range 
of  200  miles,  it  will  take  less  than  a  minute  for  the  high--  vehicles 
to  impact  and  at  most  several  minutes  for  the  low-/'  vehicles.  There¬ 
fore,  the  gravity  gradient  due  to  the  oblatencss  oi  the  eatth  is  negli¬ 
gible,  and  the  term  D  can  be  omitted  in  the  investigation  of  endoatmos- 
pheric  trajectory  prediction.  The  effect  of  the  earth's  rotation  rate 
tc  is  also  negligible,  except  for  small  deviations  that  are  observed 
during  the  last  10,000  ft  before  impact.  The  impact  point  is  defined 
in  this  report  as  the  point  at  which  a  trajectory  reaches  an  altitude 
of  10,000  ft. 

Theoretical  considerations  and  numerical  results  obtained  indicate 
that  the  term  £  in  Eq.  (1)  is  negligible  for  endoatmosphcre  prediction, 
and  this  is  especially  true  for  high-/  missiles.  An  approximate  differ¬ 
ential  equation  describing  the  ballistic  trajectory  takes  the  form 

X  »  AX  ♦  g  (2) 

Moreover,  if  as  can  be  considered  ss  zero,  then  the  differential  equation 
is  simplified  further  and  becomes 
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The  above  equation  generates  predicted  trajectories,  which  are  very  close 
to  the  actual  trajectories  down  to  50,000- ft  altitude  from  an  altitude  of 
300,000  it.  Hence,  if  the  intercept,  altitude  is  higher  than  50,000  ft, 
the  differential  F.q.  (3)  is  a  good  approximation  to  the  equation  of  motion 
for  both  high-/3  and  low-/i  missiles.  For  prediction  of  the  trajectory  down 
to  10, OOti -ft  altitude,  the  above  differential  equation  is  still  a  good 
mathematical  model  for  bigh-/3  missiles. 

The  density  of  the  atmosphere  changes  in  a  complex  manner;  an  expo¬ 
nential  curve  was  used  to  approximate  the  dens ity- al t i tude  curve.  As 
Fig.  2  shows,  this  curve  does  not  match  exactly  with  the  U.S.  Standard 
Atmosphere  1962 However,  for  the  prediction  of  an  impact  point,  this 
exponential  model  is  sufficiently  accurate. 

Several  characteristics  of  Eq.  (1)  in  endoatmosphere  are  discussed 
in  the  following  sections,  and  the  sensitivity  of  impact  points  to  initial 
values  is  mentioned  in  Sec.  V. 

A .  Characteristics  of  th  Ball  i  stic  Trajectories 
in  Endoatmosphere 

As  shown  by  the  numerical  results  (  see  Fig.  3),  pioj  ectioris  of  the 
trajectories  on  the  x-y  plane  are  almost  straight  lines.  If  the  initial 
conditions  are  the  same  for  trajectories  with  different  constant  values 
of  /?,  then  the  x-y  projections  of  their  trajectories  will  lie  on  top  of 
each  other  with  the  hi  gh-/?  mi  ssi  ies  flying  further  than  the  low-/?  missiles 
The  impact  points  lie  on  a  straight  line  in  the  x-y  plane  regardless  of 
/?  values. 

From  Fig.  3  it  can  be  seen  that  the  projections  of  ballistic  trajec¬ 
tories  on  the  x-v  plane  do  not  differ  very  much  for  different  values  of 
However,  the  projections  of  ballistic  trajectories  on  the  z- x  plane 
differ  slightly  for  different  values  of  /3  (see  Fig.  4). 

The  sensitivity  of  the  impact  point  prediction  to  the  ballistic  co¬ 
efficient  /?  is  a  function  of  the  value  of  /?;  for  low  /?  the  sensitivity 
is  large,  and  for  high  /?  the  sensitivity  becomes  small.  Therefore,  it  is 
rather  important  to  detect  whether  a  target  missile  has  low  /?  or  high  fi. 

In  th<  high-/?  case,  it  is  possible  to  predict  the  impact  point  with  high 
accuracy.  Any  inaccuracy  can  be  made  smaller  by  re-estimating  the  bal¬ 
listic  coefficient.  The  major  effect  of  /?  on  the  trajectory  occurs  at 
an  altitude  of  less  than  150,000  ft.  Since  /?  comes  into  the  differential 
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FIG.  2  DENSITY  OF  THE  ATMOSPHERE 


FIG.  3  BALLISTIC  TRAJECTORIES  ON  THE  x-y  PLANE,  CASE  1 
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equations  in  the  form  /'//?  and  p  is  a  verv  small  value  for  high  altitude, 

/  / p  is  not  a  significant  term  unless  the  altitude  is  comparatively  low 
(i.e.,  less  than  150,000  ft). 

B.  Influence  of  £  in  Endoatmosphere 

In  general,  the  shape  of  the  /3-altitude  graph  is  pa  rabo  1  i  c- 1  i  ke  and 
has  a  maximum  value.  Several  examples  are  shown  in  Figs.  5  and  6.  Ex¬ 
perimental  computations  have  been  conducted  by  using  the  minimum,  average, 
and  maximum  values  for  /?.  These  results  were  then  compared  with  the  exact 
solution.  The  predicted  trajectories  and  impact  points  for  some  repre¬ 
sentative  cases  are  shown  in  Figs.  3,  4,  7,  8,  9,  and  10.  When  the  min¬ 
imum  value  (constant)  was  used,  the  numerical  results  turned  out  to  be 
quite  different  from  the  exact  solution.  On  the  other  hand,  if  the 
maximum  value  (constant)  of  /?  is  used,  the  deviation  from  the  exact  solu¬ 
tion  is  not  too  large.  For  the  case  of  high  /?,  the  maximum  value  will 
give  the  impact  point  without  a  significant  error.  However,  the  approxi¬ 
mate  value  of  /?  will  cause  an  error  in  the  impact  time. 

If  a  target  missile  is  known  to  have  a  characteristic  of  high  /?,  it 
is  permissible  to  calculate  the  impact  point  by  using  a  predicted  maximum 
value  of  6  or  a  value  slightly  smaller  (by  10  to  20  percent)  than  the 
maximum  value  of  /?.  If  a  target  missile  has  a  characteristic  of  low  /?, 
the  prediction  of  the  impact  point  will  be  more  difficult  than  that  for 
a  high-/?  missile  (refer  to  Figs.  3  and  4). 

Consider  the  x-y  projection  of  the  trajectory.  The  impact  points 
for  high-/?  ballistic  coefficients  are  very  close  to  each  other  even 
though  the  /?  values  are  different  Some  examples  where  £  ranges  from 
1000  to  5000  lb/ ft”  show  that  the  deviations  are  5000  ft  in  the  x 
direction  and  3000  ft  in  the  y  direction  at  the  impact  point  (refer  to 
Figs.  3  and  4). 

For  low-/?  and  high-/?  missiles  having  the  same  initial  conditions 
and  impacting  on  the  surface  of  the  earth,  the  impact  points  for  the 
low-/  missiles  lie  on  the  x-y  projection  of  a  trajectory  for  a  high-/5 
missile;  in  other  words,  the  projection",  of  low-/?  missiles  on  the  x-y 
plane  arc  shorter  than  those  of  hi gh-/3  mi ssi  1  es.  Moreover,  the  x-y  pro¬ 
jections  of  these  trajectories  are  almost  straight  lines. 
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FIG.  •  BALLISTIC  TRAJECTORIES  ON  THE  z-x  PLANE,  CASE  2 


Projections  of  trajectories  on  the  z-x  plane  for  different  values 
of  ft  are  again  almost  the  some  down  to  an  altitude  of  100,000  ft.  Below 
100,000  ft,  the  z-x  projections  start  separating  and  produce  comparatively 
significant  differences  at  impact.  Some  examples  show  that,  because  of 
different  ft  values  (from  200  to  5000  lb/ ft 2),  projections  of  trajectories 
on  the  x-y  plane  differ  almost  50,000  ft  in  the  x  direction  anu  iO.OOO  ft 
in  the  y  direction  (refer  to  Figs.  3,  4,  7,  and  8). 

Naturally  the  deviations  differ  according  to  the  different  values 
of  the  initial  conditions.  If  the  initial  altitude  is  low  and  the  descent 
speed  z  is  high,  then  the  deviation  due  to  the  different  values  of  ft  is 
not  large.  For  comparison  purposes  one  example  of  this  kind  is  shown  in 
Figs.  9  and  10. 

C.  Effect  of  Nonlinearity  in  Endoatmosphere 

For  endoatmospheric  prediction,  the  nonlinear  term  C  is  negligible 
if  there  is  no  significant  change  in  values  of  r  and  V.  In  the  endoatmos¬ 
phere  case  the  maximum  deviation  of  r  wi 1 1  be 


r0  -  r( impact) 
r(  impact) 


300,000) 


where  a  is  the  radius  of  the  earth.  Hence,  the  value  [l/r|  -  l/rJ]in  the 
term  C  is  negligible  compared  to  1/r^  for  the  endoatmospheric  missile 
trajectories. 

On  the  other  hand,  the  value  (V0  -  V]  in  the  ttrra  C  is  net  necessarily 
negligible  compared  to  t*0,  Since  a  high-/5  missile  does  not  slow  down  sig¬ 
nificantly,  the  value  (Vfl  -  V]  is  negligible  compared  to  *0.  TV*  velocity 
of  a  low-/5  missile  changes  its  velocity  much  more  than  that  of  a  high-/? 
missile.  As  s  result,  the  value  "  V]  is  not  negligible,  and  for  some 
low-/?  missile* 


beresc*  0.S  or  greater. 


In  conclusion,  for  endo a tmosphe r i c  prediction  the  term  C  is  negli¬ 
gible  for  high-/?  missiles,  and  the  term  C  should  be  handled  carefully  for 
low-/?  missiles.  In  order  to  illustrate  the  effects  of  the  term  C  on 
prediction,  trajectories  of  high-/?  and  1  ow-/5  mi ssi  1  es  are  shown.  Figures 

11  and  12  and  Table  I  show  the  high-/?  case,  and  Table  II  shows  the  low~/3 
case. 

D.  Influence  of  in  Endoatmosohere 

For  endoatmospher ic  prediction,  the  influence  of  the  earth  rotation 
rate  u'  is  much  more  significant  than  that  of  the  eccentricity  e.  An  ex¬ 
ample  of  a  hi gh-,3  mi ssi 1 e  in  the  endoatmosphere  is  shown  in  Figs.  11  and 

12  and  Table  I,  With  and  without  consideration  of  o>,  the  flight  time 
difference  is  less  than  0.5  sec  and  the  deviation  of  the  impact  points  is 
about  6400  ft  (5000  ft  in  the  x  direction  and  4000  ft  in  the  y  direction). 

In  this  example,  it  takes  about  43  sec  to  impact.  A  fixed  point  on 
the  equator  moves  about  10  nautical  miles  during  these  43  seconds.  Then 
why  is  the  deviation  of  the  impact  points  with  end  without  consideration 
of  a;  about  1  nautical  mile  rather  than  10  nautical  miles?  The  answer  to 
this  q  .-scion  is  straightforward.  The  velocity  of  a  target  is  measured 
with  respect  to  the  moving  coordinate  system,  which  is  fixed  to  the 
earth  at  the  radar  site  and  rotates  with  the  earth.  Therefore,  the 
deviation  of  impact  points  with  and  without  u  is  not  caused  by  the  motion 
of  the  radar  site  but  is  mainly  caused  by  the  effect  of  the  coriolis  term 
in  Kq .  ( 1 ) . 

Another  example  is  shown  in  Table  II.  This  is  the  case  of  a  low-/3 
missile  in  the  endoatmosphere .  The  deviation  of  the  impact  points  with 
and  without  consideration  of  cc  is  about  5000  ft  (4000  ft  in  the  x  direc¬ 
tion  and  3000  ft  in  the  y  direction).  Since  this  is  a  low-/?  missile,  it 
takes  about  74  sec  to  impact,  which  is  about  80  percent  longer  than  the 
time  required  for  the  high-/?  missi  le  with  the  same  initial  conditions. 

The  low-/3  missile  takes  a  longer  time  to  impact  than  the  high-/?  missile, 
yet  the  deviation  of  impact  po-nts  for  the  iow-/$  missile  with  and  without 
consideration  of  u  is  smaller  than  that  for  the  high-/?  missile.  This  is 
because  the  main  contribution  •>f  oj  is  the  coriolis  term,  which  is  pro¬ 
portional  to  the  vector  produet  oj  *  F  of  the  earth  rotation  rate  and  the 
velocity  of  the  missile.  Since  o>  is  constant,  the  effect  due  to  coriolis 
term  depends  on  the  magnitude  of  F.  The  higher  the  velocity  of  a  missile, 
live  arger  the  deviation  that  will  occur. 
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FIG.  11  INFLUENCES  OF  •,  w,  AND  £  (Endootmosph*ro,  /3  =  /30),  x-y  PLANE 
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E.  Effect  of  Eccentricity  in  Endoatmosphere 

Eccentricity  e  comes  into  the  differential  equation  of  motion  as  the 
correction  introduced  into  the  gravitational  force  term  due  to  the  oblate¬ 
ness  of  the  earth.  The  deviations  of  impact  points  with  and  without  con* 
sideration  of  the  eccentricity  e  are  shown  in  Figs.  11  and  12  and  in 
Tables  I  and  II. 

Figures  11  and  12  and  Table  I  show  the  case  of  a  high-/?  missile  in 
endoatmosphere.  The  deviation  is  about  220  ft  (180  ft  in  the  x  direction 
and  130  ft  in  the  y  direction).  There  is  no  difference  in  impact  time. 

Table  II  shows  the  case  of  a  low-/?  missile  in  endoatmosphere.  The 
deviation  is  about  100  ft  (100  ft  in  the  x  direction  and  22  ft  in  the  y 
direction).  The  difference  in  impact  time  is  about  0.03  sec. 

In  conclusion,  the  effect  of  the  eccentricity  e  is  entirely  negligible 
for  trajectory  prediction  in  the  endoatmosphere. 


approximated 

(4) 

(5) 

in  endoatmospheric  trajectory  prediction. 

If  we  consider  p//?,  r,  and  V  to  be  piecewise  constant,  then  the 
matrix  4  and  the  term  C  become  piecewise  constant.  Hence,  it  is  possible 
to  find  a  piecewise  closed- form  solution  for  Eqs.  (4)  and  (5).  The  accu¬ 
racy  of  the  solution  depends  upon  the  integration  step  size  and  the  time 
interval  during  which  p// 3,  r,  and  V  are  kept  constant.  Experimental 
computations  were  performed  by  taking  five  time  intervals,  namely,  1,  2, 
3,  4,  and  5  sec.  The  results  are  shown  in  Tables  III  and  IV. 

The  purpose  of  obtaining  a  closed- form  solution  is  to  shorten  the 
computation  time  to  predict  the  missile  trajectory.  The  approximation 
described  above  produces  some  inaccuracies.  If  the  inaccuracies  can  be 
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tolerated,  a  closed- 1'orm  solution  should  be  used  in  order  to  reduce  the 
computation  time. 

As  a  reference,  it  may  be  helpful  to  give  an  approximate  computation 
time  to  solve  the  differential  Eq .  (1)  numerically.  If  the’  iteration  step 
size  At  is  taken  as  100  msec  and  if  one  iteration  does  not  exceed  0.2  msec 
then  the  numerical  calculation  of  a  missile  trajectory  for  a  50-see  flight 
requires  less  than  0.1  sec  of  computation  time. 


IV  EXOATMOSPHERIC  PREDICTION 


In  this  report,  the  exoatmosphere  is  defined  as  the  space  above  alti¬ 
tude  300,000  ft.  Since  the  characteristics  of  the  motion  of  ballistic 
missiles  in  the  endoatmosphere  and  exoatmosphere  are  significantly  dif¬ 
ferent  it  is  very  meaningful  to  obtain  schemes  of  predicting  trajectories 
separately . 

As  discussed  in  Sec.  Ill,  the  equation  of  motion  of  missiles  in  endo¬ 
atmosphere  is  described  by  Eqs.  (2)  or  (3).  The  equation  of  motion  of 
missiles  in  exoatmosphere  can  be  approximated  as 

i  =  AX  *  B  ♦  C  (6) 

Projections  of  the  trajectories  in  exoatmosphere  on  the  x-y  plane 
are  almost  straight  lines.  Projections  of  trajectories  on  the  z-x  plane 
have  the  shape  of  an  ellipse  or  a  parabola.  One  example  is  shown  in 
Figs.  13  and  It  and  Table  V. 

The  values  of  /. g  is  1.488  *  10” 7  lb  ft” 5  at  altitude  300,000  ft 
1/i.g  =  7.6474  y  10 ~ 2 1  b  ft'J  at  the  earth  surface).  Therefore,  the 
element  (.». g/ 2/5)1’  has  very  little  influence  on  the  solutions  regardless  of 
the  value  of  /?.  The  classification  of  missiles  is  no  longer  meaningful 
in  exoatmospher  ic  prediction  problems.  I  f  we  consider  the  term (£g/2/S)l' 
to  be  negligible,  then  bq.  (6)  is  simplified  as 

X  -  AX  ♦  ft  ♦  £  ,  (7) 

whe»  e 


n 


FIG.  13  INFLUENCES  OF  •,  «,  AND  C  (Exoo»roo*ph*r«),  x~y  PLANE 


FIG.  14  INFLUENCES  OF  •,  <•>,  AND  £  (Fxoatmo •(*«•),  z-x  PLANE 
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In  the  following  sections,  we  will  discuss  the  degree  of  the  accuracy 
of  the  above  approximations  in  Eqs.  (6)  and  (7). 


A.  Effect  of  Non  1 i nea ri ty  . i n  Exoatmosphere 

As  discussed  before,  the  nonlinear  term  .C  i s  negl  igible  for  high-/3 
missiles  in  endoatmosplie re.  However,-  the  term  C  should  be  handled  care- 
fully  for  low-/?  missiles  in  endoatmosphere% 

The  effect  of  nonlinear  term  C  (or  C)  is  not  negligible  for  any  mis¬ 
sile  in  exoatmosphere.  This  is  due  to  the  significant  changes  in  the 
values  of  r  and  V'  during  the  flight.  One  example  is  shown  in  Figs.  13 
and  14.  If  the  flight  time  is  reasonably  small,  then  tfyc  changes  in  the 
values  of  r  and  V  are  negligible  even  in  exoatmosphere.  Hence,  one  way 
of  handling  the  term  C  (or  Q)  is  to  make  it  piecewise  constant.  This 


approximation  allows  us  to  reduce  Eqs.  (b)  and  (7)  to  linear  differential 
equations  with  constant  coefficients.  As  a  result,  it  is  possible  t«.  ob¬ 
tain  a  piecewise  closed-form  solution.  This  approach  is  discussed  in  the 
next  section. 

B.  Approximation  of  Nonlinear  Term  in  Exoatmosphere 

In  the  Nike-£eus  system,  prediction  of  the  trajectory  of  a  reentry 
vehicle  is  based  on  an  analytical  closed- form  solution  of  an  approximate 
set  of  equations  of  motion.  4  In  order  to  obtain  an  analytic  solution  to 
the  equations,  the  effect  of  gravity  is  omitted,  and  the  resulting  predic¬ 
tions  are  corrected  for  gravity.  This  is  one  way  of  approximating  the 
original  differential  equation  to  find  a  closed- form  solution. 

Another  approach  is  to  divide  the  total  flight  time  into  several 
intervals  and  to  find  a  closed-form  solution  for  each  interval.  In  other 
words,  term  C  is  eliminated  and  the  initial  values  of  r  and  V'  [only  r  in 
Eq.  (7)]  arc  used  for  a  certain  time  interval  [0,  At].  The  values  of  r 
and  V  are  recalculated  by  using  the  state  values  at  time  At,  and  these 
revised  constant  values  of  r  and  V  are  used  for  ca'culating  the  trajectory 
for  the  next  time  interval  [At,  2At].  This  same  procedure  is  continued 
until  impact  is  reached. 

This  idea  is  demonstrated  for  an  exoaimospheric  trajectory,  and  the 
results  are  tabulated  in  Table  V.  The  deviations  of  the  predicted  posi¬ 
tion  from  the  exact  value  after  240-sec  flight  are  the  following: 


50,000  ft  - .  without  term  C 

20,000  ft  .  60-sec  correction  of  r  and  V  in  C 

10,000  ft  - .  30-sec  correction  of  r  and  V  in  C 

5,000  ft  -  10-sec  correction  of  r  and  V  in  C 


The  more  frequently  corrections  are  made,  the  more  accurately  we  will 
obtain  solutions.  The  suitable  number  of  intervals  for  dividing  the  total 
time  depends  on  the  error  constraints. 

If  the  total  time  T  is  divided  into  N  intervals  of  At,  then  the 

nAf  <  t  <  (n  +  1  )Af 


original  equation  is  approximated  as 


AX  At)  *  B 

ri — n  —n 


for 


0,  1 . — 

At 


3? 


and  till'  initial  conditions  are  defined  as 

*o(0)  h  . 

T 

Xn(n&t)  =  Xn_1[(n  +  l)At],  n  =  1,  2,  ....  — 

In  the  case  of  exoatmospheric  trajectories,  /?  is  negligible;  hence, 

An  and  Bn  can  be  considered  to  be  constant. 

C.  Influence  of  u.'  in  Exoatmosphere 

An  example  of  a  missile  trajectory  in  excatmosphere  is  shown  in 
Figs.  13,  14.  and  Table  V.  The  deviation  of  the  impact  points  (compared 
at  the  same  time  rather  than  at  the  same  altitude)  with  and  without  con* 
sideration  of  a’  is  about  75,000  ft  (8,000  in  the  x  direction,  25,000  ft 
in  the  y  direction  and  70,000  ft  in  the  z  direction)  after  240  sec  of 
f 1 i ght. 

The  velocity  of  the  missile  in  this  example  is  much  larger  than  that 
in  endoatmosphere  examples  shown  previously.  Hence, the  effect  of  the 
coriolis  term  is  greater. 

Although  the  effect  of  a  is  negligible  for  endoatniospheric  trajectory 
predictions,  the  effect  is  very  significant  for  exoatmospheric  trajectory 
predictions. 

D.  Influence  of  Eccentricity  in  F.xoaimosphere 

Figures  13,  14,  and  Table  V  show  the  case  of  a  missile  in  exoatmos¬ 
phere.  The  deviation  of  impact  points  with  and  without  consideration  of 
the  eccentricity  c  is  about  2100  ft  (200  ft  in  the  .r  direction,  500  ft  in 
the  y  direction  and  2000  ft  in  the  z  direction).  The  eccentricity  e  is 
negligible  (with  certain  reservation)  for  the  cases  in  exoatmosphere. 

The  main  objective  of  neglecting  the  term  D  containing  the  eccentric¬ 
ity  e  is  to  simplify  the  differential  equation  in  order  to  obtain  a  closed- 
form  solution.  For  exoatmospheric  trajectories,  the  deviation  can  be  as 
large  as  0.5  nautical  miles  after  a  240-sec  flight.  If  the  tolerance  of 
the  error  is  several  miles,  then  the  eccentricity  e  can  be  considered 
as  zero.  If  eccentricity  e  cannot  be  considered  as  zero,  then  it  is 
clearly  difficult  to  find  a  closed-form  solution. 
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One  way  to  overcome  this  difficulty  is  to  find  an  efficient  numerical 
integration  technique.  It  is  feasible  to  obtain  a  solution  of  the  dif¬ 
ferential  equation  (1)  by  using  about  0.5  sec  of  computer  time  on  a 
present-day  computer  (e.g.,  UN1VAC  1108).  The  flight  time  for  these  exo- 
atmospheric  cases  is  on  the  order  of  5  minutes  or  more;  therefore,  0.5  sec 
can  be  well  justified  for  the  computer  calculations. 

It  is  also  possible  to  neglect  the  effect  of  eccentricity  e  and  to 
simplify  the  differential  equation  so  that  the  closed-form  solution  can  be 
found.  In  order  to  support  the  above  statement,  it  is  usefc’  to  tabulate 
the  state  values  at  60  sec  after  the  initial  time  for  trajectories  with 
and  without  consideration  of  the  eccentricity  e.  The  last  two  rows  of 
Table  V  show  that  the  deviation  is  about  160  ft  after  60  sec. 
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V  SENSITIVITY  OF  IMPACT  POINTS  TO  INITIAL  VALUES 
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The  examples  considered  previously  are  used  again  for  the  sensi¬ 
tivity  analysis.  In  both  endoatmosphere  and  exoatinospiie re  cases,  10  per¬ 
cent  and  20  percent  errors  are  independently  introduced  into  each  initial 
value.  The  propagation  of  each  error  to  the  final  values  is  evaluated 
by  integrating  the  differential  Eq.  (1)  numerically  and  b\  using  the 
relationship  in  Eq.  (8).  These  results  are  shown  in  Tables  VI  through 
.X.  The  values  in  parentheses  in  these  tables  are  theoretical  results 
by  using  the  relationship  in  Eq.  (8).  According  to  Bef.  2.  the  estima¬ 
tion  errors  will  become  about  2  percent  after  5-sec  of  filtering.  Hence, 
this  sensitivity  analysis  will  give  better  results  for  more  realistic 
cases. 


This  coarse  sensitivity  analysis  gives  a  good  indication  ol  the  propa¬ 
gation  ol  errors  in  the  initial  values.  In  the  example  of  exoat.mosphere, 
the  sensitivity  analysis  and  the  numerical  integration  agree  very  well,  [n 
the  case  of  endoatmosphere,  the  analysis  and  the  numerical  integration  match 
very  well  in  most  eases,  but  the  cases  where  errors  exist  in  :(0>  and  :(0) 
do  have  significant  deviations.  Errors  in  v(0),  y(0),  .itO),  and  v(0)  propa¬ 
gate  in  the  manner  expressed  in  Eq.  (9).  Hie  position  errors  propagate 
without  any  amplification  and  have  little  influence  on  the  velocity.  The 
velocity  errors  propagate  without  any  amplification  on  t  lie  velocity  itself, 
hut  they  have  a  signilirant  effect  on  the  position  errors. 
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VI  CONCLUSION 


Hie  differential  Eq.  (1)  is  a  good  mathematical  model  of  the  bal¬ 
listic  motion.  Without  any  approximation,  it  seems  hopeless  to  find  a 
closed- form  solution  of  Eq,  (1).  The  only  way  to  find  a  solution  of 
Eq  (l)  is  numerical  integration.  As  a  result,  it  requires  a  signifi¬ 
cant  amount  of  computation  time.  This  memorandum  describes  a  simplifica¬ 
tion  of  the  differential  Eq.  (1)  of  the  ballistic  trajectories.  The  main 
purpose  of  an  approximation  is  to  obtain  a  closed-form  solution. 

The  problems  are  divided  into  two  domains,  namely  the  endoatmospheric 
problem  and  the  exoatmosphe ri c  problem.  The  endoatmospheric  problem  is 
again  divided  into  two:  namely,  the  high-/?  case  and  the  low-/?  rase.  In 
each  case,  the  influences  of  the  eccentricity  e.  the  rotation  rate  &>,  the 
ballistic  coefficient  2-  and  the  nonlinearities  are  considered.  A  summary 
of  the  influences  is  shown  in  Table  X.  In  exoatmospheric  prediction  prob¬ 
lems.  the  earth  rotation  rate  oj  and  the  nonlinear  te'  i  C  should  be  treated 
carefully,  and  in  endoatmospheric  prediction  problem,  the  ballistic  co¬ 
efficient  2  should  be  handled  properly.  The  effect  of  the  ballistic 
coefficient  is  very  significant  on  the  trajectory  at  low  altitudes  (e.g., 
for  impact  point  and  impact  time  prediction).  Further  research  effort 
should  be  oriented  toward  improving  the  estimation  and  prediction  of 
ballistic  coefficients. 

Future  work  on  the  prediction  problem  is  to  obtain  closed- form 
solutions  of  Eq*  (3)  and  (?)  Onr  possible  ws>  is  to  find  a  piecewise 
closed- form  solution  over  a  suitable  time  interval  by  taking  constant 

*V 

values  of  2  and  C  *n  Eq*  (3)  and  (7).  respectively. 


Table  X 

MM\s\tn  of  i\h.i or  c.  *.  c.  and  ,1 
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APPENDIX 


The  derivation  of  the  differential  F.q.  (1)  is  discussed  in  this 
appendix. 

If  Kap  is  the  absolute  ar  .deration  of  a  missile  P  (it  is  considered 
to  be  a  particle)  in  a  reference  frame  IV  and  a  is  the  mass  of  P,  then 
the  inertia  force  F  acting  on  P  in  N  satisfies 


(A.  1) 


The  reference  frame  N  is  a  reference  frame  in  which  the  center  C  of 
the  earth  and  the  earth's  axis,  line  NS.  are  fixed  su< h  that  the  angular 
velocity  of  the  earth  E  with  respect  to  the  reference  N.  ' .  is  given  by 


cu  ■  un 


where  a  •  2*  rad/day  and  n  is  a  unit  vector  parallel  to  line  NS  This 
reference  frame  N  is  a  good  approximation  to  a  Newtonian  reierence  frame. 
From  Fig.  I , 

\  A  A 
t  »  xi  ♦  yj  *  tk 

AAA 

where  i,  j  and  A  are  unit  vectors  defined  in  Fig.  1  and  t.  v,  and  i  are 

f.  A  \ 

the  measure  numbers  in  the  direction^  i,  snd  It.  respectively,  The 
velocity  of  P  with  respect  to  the  reference  £  is  then  expressed1  as 


.  . " 
xi  *  y;  *  tk 


The  accelerations  *<£  at  time  I*  are  related  by 

V  .  t«*  .  V*  *  2*  »  ¥ 


U.1) 


t*«a»  tk'  tin*  i‘l  with  n  t. 
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where  P *  is  the  point  fixed  in  E  that  coincides  with  P  at  time  t*.  and 
•'aP  is  \lled  the  coi  n  c  i  den  t- po  i  n  t  velocity  and  the  vector  2to  *  is 

called  the  coriolis  acceleration  of  P  for  the  reference  frames  E  and  N. 

The  coincident  point  velocity  u ap  satisfies 


l'apm  v  *0°  *a  £xr+a,x(wxr) 


*ac  = 


0  ,  ya£ 


The  acceleration  of  P  with  respect  to  £  is 


£  A  iP 


Therefore, 


A  ft 

xi  *  yj  *  ik  *  OJ  x  (co  *  _r)  ♦  2o>  * 


Since  a  and  £  are  expressed  as 


a;  =  re  sin  ,A  +  a;  cos  /x] 
ft  ft  ^ 

r  =  r,i  *  rJ  *  rik 


r  y  =  v  -  ft  sin  (^i  -  /zc) 


z  *  R  cos  (//  -  ^.c) 


The  last  two  terms  of  Eq.  (A3)  are  written  as 


(A.  3) 


(A. 4) 


x  tvp 


-oArti  ♦  or[(r  cos  fx  *  rt  sin  fx)  cos  fx  -  ry)j 

+  a;2[(r  cos  ).i  *  rt  sin  fx)  sin  fx  - 

<e(  j  cos  ^  -  y  sin  fx)  i  +  t  oj  sin  fxj  -  x  a>  cos  fxk 
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Therefore, 


A 

rM  sin  /4  cos  fi}j 
r  {  sin  /»)  sin  ^  Yk  (A.  5) 

Let  us  now  consider  the  left-1,  md  side  of  Eq.  (A.  1 ) .  The  force  £ 
acting  on  a  missile  P  is  divided  into  two  elements,  namely  the  drag 
force  F_d  and  the  gravitational  force  Ff ;  hence, 

f  «  Fd  *  Fg  (A. 6) 

The  drag  force  per  unit  mass  acting  on  a  body  is  given  approximately 
by  the  equation 


A 

N (P  ~  i'x  *  2 oil  cos  fx  ~  y  sin  /x)}  t 

♦  {y  +  2x  w  sin  n  +  o.2(rf  cos  ^x  + 

♦  {'i  -  2x  a  cos  u  *  a.!(  r  cos  u  * 

i  ^ 


.  at  ;?vpifv'p 

r.  I  L  I  L 


(A.  7) 


Next,  iet  us  find  the  expression  of  the  gravitational  force. 

The  gravitational  potential  at  P  due  to  the  earth  E  is  expressed”  as 

B‘"  ■  ; |«  -  ^7  tu  -  <r,  •  *  '„>)]■*  ••• 

where  I  is  the  moment  of  inertia  about  the  line  OP  and  I,,  I  ,  and  f 

i  1  fl 

are  the  principal  moments  of  inertia  of  the  earth  E  of  the  mass  m. 

Since  the  gravitational  force  per  unit  mass  is  the  gradient  of  the 
gravitational  potential,  the  equation 


Fg/m  =  Y/Ue/p 
holds,  where  the  operator  V  is  defined  as 


V 


1  —  +  A  JL  i  n  — 

oTj  <‘X2  YX  3 
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1'int  vectors  r  and  n  are  described  in  terms  ol'  t,  j,  and  k  as  tollows 


2 CM  (a\2  .  * 

-  ~J  (;)  s 

-  |"y  cos  /Li  +  2  sin  /x  +  ft  s  i  n  /x  c  J 


sin 


If  a  6  x  1  vector  .Y  is  defined  as 


XT  =  [x,y,  z,  i,y,  2] 

then  Eq.  (A. 10)  is  written  as  Eq.  (1)  in  the  main  text. 
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